Polarisation control of optically pumped terahertz lasers 
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Optical pumping of excited exciton states in semiconductor quantum wells is a tool for realisation 
of ultra-compact terahertz (THz) lasers based on stimulated optical transition between excited (2p) 
and ground (Is) exciton state. We show that the probability of two-photon absorption by a 2p- 
exciton is strongly dependent on the polarisation of both photons. Variation of the threshold power 
for THz lasing by a factor of 5 is predicted by switching from linear to circular pumping. We calculate 
the polarisation dependence of the THz emission and identify photon polarisation configurations for 
achieving maximum THz photon generation quantum efficiency. 

PACS numbers: 78.67.-n,78.67.De,71.35.-y,78.45.+h, 78.66. Fd, 79.20.Ws, 78.47.da 



Introduction.-¥iyicitoT\s in nanoscale semiconductor ma- 
terials exhibit low-energy excitations in the range of the 
exciton binding energy, analogous to inter-level excita- 
tions in atoms, yielding infrared and terahertz (THz) 
transitions. Thus excited exciton ladder states represent 
a natural system for generating THz radiation and co- 
herence. The demand for development of new compact 
and efficient coherent terahertz radiation sources is cur- 
rently rapidly increasing, due to ever growing range of 
very diverse technological applications in the relatively 
little-explored THz spectrum of radiation Towards 
this goal recently a new scheme of a microcavity based 
polariton triggered THz laser (THz vertical cavity surface 
emitting laser (VCSEL)) has been proposed by one of the 
authors 0, whereby the 2p dark quantum well (QW) ex- 
citon state is pumped by two-photon absorption using a 
cw laser beam. 

In this Letter we theoretically demonstrate polarisa- 
tion control of THz emission and of the quantum effi- 
ciency for THz photon generation. We consider a THz 
VCSEL proposed in Ref.( 0]), where the pump beam is 
split in two. Each of the split beams goes through a po- 
lariser, so that the two photons pumping the 2p exciton 
do not necessarily have the same polarisation. We show 
that by rotating one of the polarisers one can switch on 
and off the THz laser. 

Using crystal symmetry point group theoretical meth- 
ods [sj we calculate the polarisation dependence of the 
optical transition matrix element for two-photon exci- 
tonic absorption in GaAs/AlGaAs quantum wells, as well 
as of the intra-excitonic 2p to Is THz transition radiative 
decay rate. This enables us to calculate the polarisation 
dependence of the quantum efficiency for THz photon 
generation and thus identify maximum efficiency regimes 
of operation. The optically pumping scheme to a 2p ex- 
citon state by two photons, each of half the energy of the 
2p exciton state, is shown in Fig. [TJ 
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FIG. 1: (Color online) Schematic energy-level diagram of two- 
photon transitions to 2p exciton states in a QW. The ground 
(Is) and excited (2p dark) discrete (bound) excitonic states 
and the exciton (unbound) continuum states are shown. The 
pumping frequency, UJpump is half that of the 2p exciton state. 
ojthz is the center frequency of the emitted THz pulse; To - 
ground (Is) state exciton spontaneous emission rate. 



Two-photon p-exciton absorption- The quasi-2D 
(Q2D) exciton wave function at the F point is given for 
narrow QWs by 

*A (r„ r;.) = ^Uf (p) (2,) <i>^* (zn) Kek (r.) <k (r^) e*"'^" 

(1) 

where wo is the unit-cell volume, S is the QW area, R 
is the centre-of-mass (c.o.m.) co-ordinate, r = re — is 
the relative motion co-ordinate, r(p, z) and p = — 
is the in-plane relative motion co-ordinate, the z axis 
is taken normal to the QW layers, a, /? are subband in- 
dices and is the a-subband envelope function of the 



conduction (valence) band; (re — r^) is the envelope 
function of the 2D exciton associated with subbands a 
of the electron and (3 of the hole; A — (n, m) is the 2D 
exciton quantum number, labelling the discrete excitonic 
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states {n = 1,2,..., |m| < n); Mck,w-uk are the periodic 
parts of the Bloch wave function for conduction and va- 
lence bands, correspondingly; the exciton c.o.m. wave 
vector, k|| R:! 0, is on the order of the photon wave vector. 
Consider the case of allowed conduction-to- valence band 
dipole optical transition at the F point. In cubic crystals 
the conservation of parity upon absorption of two pho- 
tons requires the final excitonic state to have the same 
parity as the valence band, therefore the final exciton is 
in a p-state. The TPA probability is given by: 



WTPA = j-Y.\V.ffSfiE) 

if 



(2) 



where Sf is the final density of states and the momen- 
tum, p, matrix element, Vif, between the initial and final 
states is given by [5|: 



- Ei — hu]o 



(3) 

reflecting the order of absorption of the first photon with 
polarisation vector £i, energy hwi and vector potential 
Ai and the second - with polarisation vector £2, energy 
huj2 and vector potential A2. Since the TPA is a two- 
step process, one should sum over all intermediate states 
\l) with energy Ei. The first matrix element has been 
calculated by Elliott [6] in the 3D case and for the quasi- 
2D case here considered reads: 



(?|e„.pN> = V^*A (0) (c|£a.p|«) 



(4) 



where a = 1,2 is the photon polarisation vector, 
^'a(O) is the relative motion exciton wave function, given 
by Eq.(IT]), evaluated at r = 0, and the interband matrix 
element (c| £a.p \v) is given by: 



vq J 

cell 



(r)-Vu„2 (r) 



(5) 



The second matrix element entering Eq. ^ is between 
hydrogenic-type exciton states and can be written as: 



l(/|e>.p|0 = - 



d^'rvl/f* (r) (e>.p)*f (r) 



= - / dh- C/f * (p) (z.) € (Zh) (e>.p) Uf (p) (z,) <l>f (^h) 

(6) 

where m is the free electron mass and fi^ is the reduced 
exciton mass along ^-direction in the QW plane. Intro- 
ducing a special notation for the matrix element, summed 
over the intermediate states, according to: 

Is {a, 0)^ J <f r (p) ize) {zh) (e^-p) 

X E ¥^rl^^^*r (z.) iz,) <&f (z,) 

'—^ Ij\+ IjG - n^a 



where Eq is the direct interband energy gap, the total 
matrix element can be written as: 



Vh = VS^^A^A^ [{c\ li.p \v) Is (1, 2) + (c| I2.P 1^') h (2, 1)] 

(8) 



Let us define a 2D reduced Coulomb Green's fimction: 
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where E\ is the exciton hydrogenic energy as measured 
from the conduction-band edge and = —Ea + HuJa < 0. 
A closed form of the reduced Green's function for an 
unscreened exciton in 3D (N-D) has been derived in Q 
(0) and in the 2D limit of interest is given in Q: 



1 -2p/ . 
G(p,0) = — e 



In 
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7 + 3-4 



(10) 

where a*^ is the 3D exciton Bohr radius, 7 is the Euler's 
constant and — ^^J^^^ with Eb = Ry* - the exci- 
ton binding energy. Introducing cylindrical co-ordinates 
and the overlap integral of the subband envelope wave 
functions: 

I^p = (z)| (z)) ^ jdz $r {z) <i>f {z) (11) 



and momentum matrix element along z-direction: 

P^p = (z)b, \<^t [z)) = j dz^T iz) (z) 

(12) 

the general expression for the sum over intermediate 
states can be recast as: 



+P^0jd^pU^^* ip) (£>.z)G(p,0) 



(13) 

where P = and z are unit vectors. 

For VCSEL configuration and normal incidence geom- 
etry we choose e^-Lz and therefore the second term in Eq. 
([T5| vanishes. The derivative of the Green's function can 
be easily carried out using Ea. (ITU|) . We take for the exci- 
ton relative motion wave function the 21) hy drog en atom 
wave function for bound exciton states [10|, |ll| : 
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where A^^ 



(n — 1 — |7n|) 



_ ("-5) [("-1+1™!)']' 
associated Laguerre polynomials [12| . 
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and L" (x) 
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FIG. 2: (Color online) 3D surface plots of the normalised excitonic TPA against polar angles of ei and £2 in the QW plane at 
specific phase shifts 5i,(52- (a) co-linearly polarised photons 5i — 0, ±7r; S2 — 0, ±7r; (b) 1"* linear - 2"'' a'^ circularly polarised 

linearly polarised photon 5^ = f ; ^2 =0; (d) a~^- cr+ 
a~- counter-circularly polarised photons 



photon 5i = 0; 82 



(c) cr+ circularly polarised photon - 2" 



or a - 

<5i =± 



a co-circularly polarised photons 5i 



±f;52 



(f) co-left-elliptically polarised photons Si 



±f ; (e) a+- 



Introducing polar co-ordinates (p, ip) we obtain for 
TPA to 2|)-exciton states with n — 1, m = ±1: 



/2,i(«,/3) 



Ahi, 



Jpail^a) (15) 



where Eb = „- ,1 is the exciton binding energy and the 
integral, Jp,2 {K) = _ Ki43+36^n( 3)) ^^3^3 ^ 

Substituting in Eq. (|5]) the excitonic two-photon ab- 
sorption matrix element is obtained: 



AM, 



al3 



1 



Sttivott EBa*B 



J,ff.2{c\p\v) (16) 



where we have defined effective matrix element for cubic 
crystals, using the invariance of the inter band matrix el- 
ement M = (c| p\v) under crystal point symmetry group 
transformations [ijl, 01, 01, [Hj: 



Jeff,2 = 2 ('"1 ^ ^2) \Jp.2 (fcl) - Jp,2 (k2) 
+ i [1 + {£l.£2)'] \Jp,2 ih) + Jp,2 (fc2)|' 

_c! 



(17) 



^ {(£1 X £2)' - fell' + [1 + (£l.£2)'] + fcfl'} 



^here = -^i^^^SI^. 

Our pumping scheme envisages two photons each with 
half the energy of the 2p-exciton state: hui = huj2 = 
fwj = ^ and kf = k"^ = = 2Eg-E2p ' therefore the 
first term in Eq. PT)) vanishes and from Eq. ^ we get 



for the TPA probability to 2p-exciton states in [s ^m ^]: 



■E2p) 



l + (ei-£2)' 
(18) 

where Sjp' is the final 2p-exciton density of states per 
unit area for a heavy- hole exciton (cl-hhl) and the coef- 
ficient Ktpa for an infinite quantum well, is given by: 



Ktpa 



USne^AfA 



(19) 



The photon polarisation vectors £1,62 with polar angles 
<y9i,(y52 and phase shifts 5i,52 correspondingly, lie in the 
QW plane. 3D plots of the exciton TPA probability are 
shown in Fig. [2] for different polarisations of the two 
pumping photons. 

We suggest adding an external THz cavity at the VC- 
SEL output that will filter out the linear polarisation of 
the emitted THz radiation, and will thus constitute our 
reference frame, fixing the direction of our co-ordinate 
system x-axis. We shall assume that the generated THz 
mode is X-polarised. By inspection of Fig. 2 one can see 
that maximum (5-fold) increase of the two-photon ab- 
sorption rate with respect to YY polarisation is achieved 
for hnearly XX, XX, XX, XX polarised photons (Fig. 
[2ta)). The two-photon absorption rate can increase by 
a factor of 3 for linearly-circularly or circularly-linearly 
polarised photons (Fig [2jb,c)); by a factor of 2 for both 
circularly polarised (Fig. [l](d,e)), by a factor close to 5 
(but always less than the one for linear polarisation) for 
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elliptically polarised photons (Fig. [5] (f)). Our results 
show that changing polarisation from YY to XX, pass- 
ing through circularly and elliptically polarised pumping, 
one can vary the lasing threshold by a factor of 5. 

Intra- excitonic 2p — >■ Is transition probability - We cal- 
culate next the polarisation dependence of the 2p ^ Is 
photon intra-excitonic transition rate, generating THz 
emission (Fig. [T]). We are interested in the optical transi- 
tion matrix element between initial two-fold degenerate 
state '^s with 6 = {n = 2,m = ±1) = {2,p) and final 
state ^'a with A = (n = l,m = 0) = (l,s). The matrix 
element is of the second type Eq. ^ and for normal inci- 
dence geometry (i-Lz) and exciton wave functions, given 
by Eqs. (P) , Ht]) , we obtain: 



(20) 

where we have introduced polar co-ordinates and the an- 
gular dependence is given by: ^ {(fi) = cos(/9eT*'^. The 
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Fi- 



integration over p is easily performed, giving: 
nally, the 2p — Is intra-excitonic optical transition rate 
for an infinite QW is given by: 



(2) _ 



2p->ls 



27 
512' 



THz "- 



/4^5i,(i?)$2(^) (21) 



where Sis{E) is the final (Is) state density of states and 
the one-photon THz emission coefficient is given by: 



^OPA 



167r 



Lis 



(22) 



where Athz is the THz photon vector potential, ex- 
pressed in terms of the THz emission intensity, Ithz as: 

/ sl/2 

Athz — ( ^^"^a ) , where n is the refractive index 
and LOT Hz = 

The polarisation dependence of the THz emission rate 
can be inferred from the angular dependence: for linear 
(e.g. along x-axis) polarisation of the emitted THz pho- 
ton (ejl x), <i>^ (<p) = 1, for y- linear £|| y, <i>^ {if) — and 
therefore there is no THz emission, and for circularly po- 
larised THz photon, $^ ((^) — i, the corresponding THz 
emission rate is half of the one for x-linear polarisation. 

Quantum efficiency - The quantum efficiency of THz 
radiation generation can be defined as the ratio of the 
THz photon generation rate and the two-photon absorp- 
tion rate by a 2p-exciton and is proportional to the ratio 
of the squares of the oscillator strengths, G and g, of the 
2p — Is and |0) — ?► \2p) transitions 12!], which can be 
expressed in terms of the transition probability [l3| : 



67rc^E/2^5 (2) 



(23) 



where uj2p 



h 



Using Eq. p^ and Ea. (pT|) . after some algebra one can 
obtain for the normalised quantum efiiciency: 



e'^h?L,h'P (143 + 36 In (§))' (^2,, - EuY 



where the polarisation dependence is given by: 



(24) 



P(£l,£2,£) 



COS^ {ip) 



1 + (cos 1^1 COS 1^2 + COS {ipi + 5i) COS {ip2 + S2)) 



(25) 

and Lp is the polar angle of the THz photon polarisation 
vector. 

We shall assume that the 2p exciton excited by two- 
photon absorption has a lifetime, which is long enough 
that it loses any memory of the polarisation and phase of 
the excitation, so that it can emit with any polarisation. 
We shall consider emission with one particular polarisa- 
tion (either linear or circular) and all possible choices of 
polarisation of the two pumping photons. 

The quantum efficiency polarisation dependence is 
shown in Fig. [3]for different polarisation configurations of 
the two pumping photons at a given (linear) emitted THz 
photon polarisation. The plots for circular polarisation 
of the THz radiation look exactly the same but are scaled 
down by a factor of 2 (not shown), resulting in maximum 
efficiency = 0.5. In addition to the results presented in 
Fig. [3l we should note that for counter-X (X) -linearly 
polarised pumping photons 5i — (tt) ; ^2 = tt (0) max- 
imum quantum efficiency 77 = 1 is achieved for linearly 
polarised (0 = 0,7r) and rj = 0.5 for circularly polarised 
— ^) THz emission, unconditionally, for any direction 
of the linear polarisation of the two pumping photons 
in the QW plane. Furthermore, if the THz emission is 
y-linearly polarised, the quantum efficiency 77 = 0, i.e. 
no THz radiation should be emitted in this case. Fig. 
13] shows that the maximum quantum efficiency ry = 1 
could be achieved within certain regions in the plane for 
linearly polarised THz emission for all combinations of 
linear and circular polarisations of the two pumping pho- 
tons. Note that in both (circular and linear THz emis- 
sion polarisation) cases, maximum quantum efficiency is 
achieved along YY lines for co-linearly polarised photons, 
for y-polarised first (second) photon in the case of linear- 
circular (circular-linear) polarisation, or along diagonal 
lines for co- and counter- circular-circular polarisation of 
the pumping photons. We emphasise, however, that al- 
though maximum quantum efficiency could be achieved 
both by counter- and co-linearly polarised photons, the 
quantum efficiency in the former case is constant and 
does not depend on the direction of the polarisation vec- 
tors in the QW plane, while maximum quantum efficiency 
in the latter case is obtained solely for specific directions 
of the polarisation vectors in the plane (YY) . 

In order to verify these predictions experimentally, one 
can envisage pumping of a QW structure with two laser 
beams having the same frequency (equal to a half of the 
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B/Wlr' inPSP WwUP' BJ^UiP'' A/wSI' 

"^W* '"'1^^ ^-^^W' ^^S^* 

FIG. 3: (Color online) 3D surface plots of the normalised quantum efBciency of THz photon generation against polar angles 
of the pumping photons polarisation vectors in the QW plane at different phase shifts Si, S2 at linear, = 0, ±7r polarisation 
of the emitted TIfz radiation (a) co-linearly polarised photons;(b) linear-2"'' circularly polarised photon; (c) 1"* circular 
-2"'' linearly polarised; (d) a'^-a^ or a~-a~ co-circularly polarised photons; (e) a^-a~ or a~~a^ counter-circularly polarised 
photons. 



2p-exciton resonance frequency) but different polarisa- 
tion. These two beams may be generated by the same 
laser but should propagate through different polarisers 
before focussing on the sample. In addition we suggest 
including a delay line between the two parts of the pump- 
ing beam, which would provide the phase difference of tt 
between them to obtain counter-linearly polarised beams 
for which unconditional maximum efficiency is predicted. 
The intensity and polarisation of the THz light emitted 
by the structure could be measured as a function of in- 
tensities and polarisations of the two pumping beams. 
As reference experiments one can measure the intensity 
of THz emission with one of the pump beams switched 
off. Analysing the results of such experiments one should 
bear in mind that the two photons used to generate a 
2p-exciton may originate from the same beam as well as 
from different beams. Comparing the spectra obtained 
with both beams switched on with those obtained with 
only the first or only the second beam switched on, one 
can extract the signal generated by absorption of the two 
photons coming from different beams and thus having 
different polarisations. 

Conclusions.- We have developed a theory of the two- 
photon absorption to p-exciton states in QWs and cal- 
culated the polarisation dependence of two-photon tran- 
sition probability, using crystal symmetry point group 
methods. We show that the two-photon transition rate is 
strongly dependent on the polarisation of both photons 
and our model predicts variation of the lasing thresh- 
old by a factor of 5 by switching from YY to co-linearly 
XX-polarised pumping. We calculated the polarisation 
dependence of the intra-excitonic THz emission and the 
quantum efficiency for THz photon generation. Maxi- 
mum quantum efficiency is predicted for counter-linearly 
polarised pumping photons and linearly polarised THz 
emission. Conditions for achieving maximum quantum 
efficiency for different polarisations of the pumping pho- 
tons are identified, thereby opening routes for polarisa- 
tion control of the THz VCSEL and a range of new ap- 
plications entailed from it. 
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